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ABSTRACT 

We present XMM-Newton observations of the Seyfert 1 AGN H 0557-385. We 
have conducted a study into the warm absorber present in this source, and using 
high-resolution RGS data we find that the absorption can be characterised by two 
phases: a phase with log ionisation parameter £ of 0.50 (where £ is in units of ergs cm 
s _1 ) and a column of 0.2xl0 21 cm~ 2 , and a phase with log £ of 1.62 and a column 
of 1.3xl0 22 cm~ 2 . An iron Ka line is detected. Neutral absorption is also present 
in the source, and we discuss possible origins for this. On the assumption that the 
ionised absorbers originate as an outflow from the inner edge of the torus, we use a 
new method for finding the volume filling factor. Both phases of H 0557— 385 have 
small volume filling factors 1%). We also derive the volume filling factors for a 
sample of 23 AGN using this assumption and for the absorbers with log £ > 0.7 we 
find reasonable agreement with the filling factors obtained through the alternative 
method of equating the momentum flow of the absorbers to the momentum loss of 
the radiation field. By comparing the filling factors obtained by the two methods, we 
infer that some absorbers with log £ < 0.7 occur at significantly larger distances from 
the nucleus than the inner edge of the torus. 

Key words: Galaxies: active - galaxies: Seyfert - galaxies: individual (H 0557—385) 
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1 INTRODUCTION 

Warm absorbers are clouds of ionised gas intrinsic to AGN. 
They cause absorption at soft X-ray wavelengths, and 
nearby bright Seyferts are ideal for spectroscopic studies of 
this phenomenon. Warm absorbers are a useful probe of the 
inner environments of AGN, e.g. elemental abundances, out- 
flow speeds, and degrees of ionisation. 

The importance of warm absorbers was underlined fol- 
lowing a study by Reynolds (1997) which indicated that they 
are present in about half of all type 1 AGN. The advent of 
XMM-Newton and Chandra has shown the details of warm 
absorbers at their highest resolution yet, revealing that their 
primary spectroscopic signature is a series of narrow absorp- 
tion lines (Kaastra et al 2000). 

Many warm absorbers in type 1 AGN have now 
been observed with high resolution X-ray instruments, e.g. 
NGC 5548 (Kaastra et al 2000); IRAS 13349+2438 (Sako et 
al 2001); NGC 3783 (Kaspi et al 2002); IC 4329A (Steen- 
brugge et al 2005). Observations of warm absorbers typically 
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show that there are at least two phases present, with dis- 
tinctly different ionisation parameters, although in theory 
a continuum of ionisation parameters could be present be- 
tween the observed phases. An unresolved transition array 
(UTA) of M-shell iron, first discovered in IRAS 13349+2438 
(Sako et al 2001), has been found to be an important spec- 
troscopic diagnostic tool (Behar et al 2001). 

Two main theories exist as to where warm absorbers 
originate: either as evaporation from the obscuring molecu- 
lar torus (e.g. Krolik & Kriss 2001) invoked in Seyfert unifi- 
cation models (Antonucci 1993), or as an accretion disc wind 
(e.g. Elvis 2000). A detailed study of 23 AGN by Blustin et 
al (2005; hereafter B05) indicates that the warm absorbers 
of nearby Seyferts are most likely to originate as outflows 
from the molecular torus. 

A recent hypothesis for the geometry of warm absorbers 
is that they have the form of an 'ionisation cone': this sce- 
nario involves a cone of plasma, irradiated by the nuclear 
continuum, with a shape defined by the torus. From obser- 
vations of the Seyfert 2 AGN NGC 1068, this scenario is 
suggested to account for both absorption features in Seyfert 
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1 AGN, and emission features in Seyfert 2 AGN (Kinkhab- 
wala et al 2002; Brinkmann et al 2002; Ogle et al 2003). 

H 0557-385 (also known as 3A 0557-383, EXO 055620- 
3820.2 and CTS B31.01) is a Seyfert 1 AGN at redshift 0.034 
(Fairall, McHardy & Pye 1982). The Galactic neutral col- 
umn density towards it is 4xl0 20 cm -2 (Dickey & Lockman 
1990). We report the first observations of this source with 
XMM-Newton. H 0557-385 has also been observed with 
ASCA (Turner et al 1996) and BeppoSAX (Quadrelli et al 
2003). Both these previous observations detected a warm 
absorber. In this paper we explore the details of the warm 
absorption in H 0557— 385, using both EPIC and RGS data. 



2 OBSERVATIONS AND DATA REDUCTION 

H 0557— 385 was observed with XMM-Newton (Jansen et 
al 2001) on 3rd April 2002 and 17th September 2002; the 
observation details are summarised in Table Q 

The EPIC pn data were taken in large window mode, 
and the MOS data in small window mode, both using the 
thin filter. We excluded periods of high background, which 
we identified in 5-10 keV lightcurves for the whole field of 
view outside the source region. The data were processed 
using the Science Analysis System (SAS) Version 6; the pn 
spectra were accumulated using single and double events, 
corresponding to PATTERN values of 0-4, and the MOS 
spectra were constructed using all valid events (PATTERN 
= 0-12). The count rates are below the thresholds where 
pile-up has to be considered (12 counts s -1 for the pn large 
window mode, and 5 counts s _1 for the MOS small window 
mode). Events next to bad pixels and next to the edges of 
CCDs were excluded (FLAG = in SAS). 

The EPIC source spectra were constructed by selecting 
counts within a circle of radius 25 arcsec, and the back- 
ground spectra were extracted from nearby source-free re- 
gions in circles of radius 75 arcsec. 

Response matrix files (RMF) and auxiliary response 
files (ARF) for each instrument were generated using the 
SAS. The count rates for each observation are similar, as 
shown in Table The EPIC spectra from both observa- 
tions were coadded using the method of Page, Davis & Salvi 
(2003) and binned into 45 eV channels. 

The RGS spectra were extracted using rgsproc in SAS 
V.6. The first and second order spectra and response matri- 
ces from both RGSs and both observations were resampled 
to the channels of the first order RGS1 spectrum from April 
2002. The spectra were then coadded, and the response ma- 
trices were combined, to produce a single spectrum and a 
single response matrix. The spectrum was then binned by a 
factor of 9 to improve signal to noise, giving a binsize of 0.1 
A at 20 A. 

The data were analysed using SPEX 2.00 (Kaastra et 
al 2002). In this paper we adopt the values Ho = 70.0 km 
s -1 Mpc~\ n m = 0.3 and Q A = 0.7. 



3 RESULTS 

3.1 Fits to the EPIC data 

After initial inspection, the EPIC data were fitted in SPEX 
with a number of spectral models: the fit parameters are 



listed in Table [5] and data and fitted models are shown in 
Fig. Absorption by the Galactic column density is in- 
cluded in all the fits, using the model hot, in SPEX. This 
is for a gas in collisional ionisation equilibrium, set at the 
minimum temperature (5x 10 -4 keV), and it produces a very 
similar spectrum to that of a pure photoelectric model (e.g. 
Morrison & McCammon 1983) , but includes absorption lines 
as well as absorption edges. 

We started by fitting a simple power law in the range 
2.3-10 keV (Model A) to avoid the regions where the low- 
energy absorption affects the spectrum. In Fig. 0A) the 
power law is plotted back to 0.3 keV to show the shape and 
extent of the absorption. The plot is suggestive of the pres- 
ence of a warm absorber, with the deepest trough occuring 
at ~ 0.7 keV; at this energy this will be a warm, rather than 
a neutral absorber. 

To fit the whole spectrum, a warm absorber compo- 
nent (xabs in SPEX) was added to the model (Model B). 
The xabs component models absorption by a cloud of pho- 
toionised gas, at a given column and ionisation parameter, 
where the ionisation parameter is defined as £ = L/nr 2 in 
erg cm s _1 . Here L is the 1-1000 Rydberg ionising luminos- 
ity (erg s _1 ), n the gas density (cm -3 ), and r the distance 
of the ionising source from the absorbing gas, in cm. We set 
the turbulent velocity of the gas to 100 km s - , and the 
elemental abundances to solar values. Only the column and 
ionisation parameter were left free in the fitting. 

Including the warm absorber in the model brings it to 
resemble the data far more closely, as shown in the plot for 
Model B. However, the xVdof is 946/211, and the model 
is still not a good match to the data at some energies. The 
single xabs component has created an absorption feature in 
the model which is deeper than that in the data at ~ 0.75 
keV, but the model has insufficent absorption at 1 keV, as 
examination of the data/model ratio in Model B shows. 

Another xabs component with a much higher ionisation 
parameter was added (Model C): this reduces the ionisation 
parameter and column of the first xabs component and pro- 
vides a better fit to the data at ~ 0.75 keV and 1 keV, 
as shown in Fig. 0G)- The goodness of fit is significantly 
improved, but the model still appears to require some ab- 
sorption at ~ 0.3-0.4 keV to match the data. Therefore a 
component of neutral gas (again, using hot in SPEX at the 
minimum temperature of 5xl0 -4 keV) was added at the 
same redshift as H 0557—385. The additional column den- 
sity of neutral gas, 4x 10 20 cm -2 , provides a much better fit 
to the data at lower energies, as shown in Fig. Q^D). With 
the addition of neutral gas, the column of the high-£ ab- 
sorber increases and the model fits the data more closely at 
~ 0.3-0.4 keV and ~ 0.5-1 keV. 

Now the soft X-ray spectrum is fitted well, we turn to 
the hard X-ray spectrum. An iron Ka line is visible at ~ 6 
keV and we first modelled it with a gaussian. This improved 
the x 2 by 49 for 3 extra parameters (Model E in Fig. 0. 
The energy of the line is 6.41±0.05 keV, with a FWHM of 
3301 

130 e V and equivalent width of 85^24 eV. As this line 
may originate by fluorescence in reflecting material, we next 
tried a reflection model (refJ in SPEX), which includes an 
iron Ka line, in combination with Model D. As the energy 
of the line indicates it comes from cold material, we assume 
the reflecting material is cold. The best-fit reflected fraction 
is 0.39^q' 1 o- The addition of this component improves the 
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Table 1. Observation details. The RGS count rates are for the co-added first and second orders of RGS1 and RGS2, for each observation. 
The Effective Exposure Time is the exposure time with periods of high background excluded. 



Date Observed 


Effective Exposure Time 
(ks) 


Background-subtracted Count Rate 
(Counts s — 1 ) 


03.04.2002 
17.09.2002 


pn: 8, MOS: 3, RGS: 12 
pn: 7, MOS: 9, RGS: 8 


pn: 10.29±0.05, MOS1: 3.30±0.04, MOS2: 3.32±0.04, RGS:0.125±3 X 10" 3 
pn: 10.55±0.04, MOS1: 3.63±0.02, MOS2: 3.68±0.02, RGS:0.143±3 X 10~ 3 



X 2 by 44 from Model D. This is shown in Fig. [I] as Model 
F. As the reflected fraction is small, and the addition of this 
component produces a slightly worse \ 2 than Model E, we 
conclude that the inclusion of reflection is not justified and 
that the iron Ka line alone is a sufficient addition to Model 
D. We then take Model E as the best-fit to the EPIC data. 

The 2-10 keV flux of model E is 4.3±0.1 x 1CT 11 erg s" 1 
cm -2 , compared to2xl0 -11 ergs -1 cm -2 at the time of the 
ASCA observation of Turner et al (1996) and 4xl0~ n erg 
s _1 cm -2 at the time of the BeppoSax observation of Q03. 
The 1-1000 Rydberg luminosity of model E is 4.3 xlO 44 erg 
s- 1 . 



3.2 Fits to the RGS data 

We now investigate the soft X-ray absorption features using 
the RGS, which has much higher energy resolution (E/AE 
~ 200-800) than EPIC (E/AE ~ 20-50). 

We plotted the EPIC best-fit, model E, on to the RGS 
spectra after convolving it with the RGS response; the re- 
sult is shown in Fig. H RGS and EPIC cross-calibration 
shows discrepancies of up to 15 per cent in the band 6-18 
A and up to 10 per cent in ~ 18-36 A (Fig 11 of Kirsch 
et al 2004). The spectral shape and flux of H 0557-385 are 
generally well reproduced, but this model does not match 
the RGS data completely. It overestimates the data over ~ 
8-15 A and ~ 23-33 A and is too absorbed at ~ 18-19 A. 
In order to obtain a better fit to the RGS data, we kept the 
continuum from EPIC model E fixed, and the columns and 
ionisation parameters of the absorbers in the model were left 
to vary. We consider only the 8-30 A spectrum, where RGS 
calibration and statistical quality of the data are best. The 
turbulent velocity of the gas was again set to 100 km s . 
The absorbers were set at the recession velocity of the AGN. 
The resultant best-fit proved to be in much better agreement 
with the RGS data, and is shown in Fig. [3] along with the 
data/model ratio. The main ions from each absorbing phase 
are indicated. The parameters from the EPIC model E and 
RGS fits are compared in Table. [3] 

As the RGS has better energy resolution than EPIC, 
the best-fit model will be closer to the real characteristics of 
the source. As shown in Table. [31 the parameters of the warm 
absorbers required by the RGS data are different to those 
from the EPIC data. With respect to the EPIC data, for the 
RGS data the column of the low-£ phase is lower whereas the 
column of the high-£ phase has not changed. The log £ values 
of the low-£ phase are consistent, but the log £ of the high- 
£ phase is smaller than that derived from the EPIC data. 
A much larger column of neutral gas is also required. The 
discrepancy between the EPIC and RGS column densities 
for the neutral and low-£ absorbers could be related to the 
larger uncertainty in the pn response at the lowest energies. 



Compared to Fig.|5] the model in Fig. is now a much 
better fit to the data at 8-15 A and from 17 A onwards. The 
X 2 /dof of the model shown in Fig.JSJis 413/236, whereas the 
X 2 /dof of the RGS best-fit model (Fig. ^ is 246/231. The 
discrepancy between the data and the model at ~ 18-19 A 
in Fig.|5|has disappeared. 

Fig-Elshows the models corresponding to the high-£ and 
low-f phases. The deep and broad absorption feature at ~ 
16-18 A, clearly visible in the RGS spectrum, is the signature 
of an unresolved transition array (UTA) of M-shell iron. The 
log £ values of the two phases are such that a UTA is present 
in both of them. There is much more absorption at shorter 
wavelengths for the high-£ phase than for the low-£ phase, 
as expected. 

To check the self-consistency of the RGS-derived ab- 
sorption parameters with the continuum from the EPIC fit, 
we re-fitted the EPIC spectrum with the absorption parame- 
ters fixed at the RGS values. The best fit continuum changes 
from r=1.90 to r=1.96, consistent at 3 sigma; the difference 
is well within the absolute uncertainty in the EPIC calibra- 
tion, and is much less than the typical difference resulting 
from the EPIC-RGS cross calibration. 



4 DISCUSSION 

The EPIC data tell us broadly the characteristics of the 
warm absorbers; model E tells us that the warm absorption 
is caused by two phases with different ionisation parame- 
ters. A column of neutral gas is also required to explain the 
spectrum at ~ 0.3-0.4 keV. The RGS data, at higher res- 
olution, require a smaller log £ value for the high-£ phase 
and a smaller column for the low-£ phase. A larger column 
of neutral gas is also required for the RGS spectrum. 

The Fe Ka line, with an energy of 6.41±0.05 keV, orig- 
inates in material that could be neutral, or ionised up to 
ionisation of Fe XVIII, which is consistent with this energy 
(Kallman et al 2004). This wide range of ionisation states 
corresponds to log £ < 2.3, so we cannot constrain well 
the ionisation state of the line-emitting material. Assum- 
ing it is a single line, it is resolved in EPIC with a FWHM 
of 16,000±|ggg km s" 1 , whereas the ASCA observation of 
Turner et al (1996) yielded a FWHM of 55,000+2^000 km 
s _1 . However, the width of the line could be due to con- 
tributions from multiple ionic species of iron. We find the 
equivalent width of the line to be 85+24 e V, compared with 
~ 300 eV found by Turner et al (1996). 

From optical spectroscopy, Rodriguez-Ardila et al 
(2000) distinguish 3 Ha components in H 0557-385. The 
FWHM of the narrowest component is 1035 km s , that of 
the intermediate component 2772 km s _1 and for the broad- 
est component it is 11,000 km s _1 . The last value suggests 
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Table 2. : Parameters of fits to the combined EPIC data. The Galactic column (4xf0 20 cm -2 ) is included in the models. PL norm is 
the normalisation of the power law in units of 10 52 ph s _1 kcV -1 . £ is the ionisation parameter in units of erg cm s~ x . The reflection 
fraction is defined as Q/2n, where fi is the solid angle subtended by the reflector at the X-ray source. Errors are at 90 % confidence for 
one interesting parameter (Ax 2 = 2.71). 



Model 


r 


PL norm 
at 1 kcV 


Warm absorber 
N H (10 22 cm- 2 ) 


log 5 of 
warm absorber 


Neutral N H 
(10 22 cm- 2 ) 


Fe Ko 
energy (kcV) 


Fc Ka 
a (eV) 


Fe Ka 
EW (cV) 


Reflection 
fraction 


x 2 /d.o.f 


PL X xabs 
(Model B) 


1.79±0.01 


6.91±0.13 


0.76±0.02 


0.41±0.02 












946/211 


PL X xabs X xabs 
(Model C) 


1.83±0.01 


7.48±0.20 


0.62±0.02 
0.7B±0.10 


0.07±0.0B 
2.2B±0.06 












637/209 


PL X xabs X xabs 
Xncutral N H (Modcl D) 


1.88±0.01 


8.11±0.20 


0.60±0.02 
1.11±0.20 


0.37±0.06 
2.30±0.07 


0.04±0.01 










487/208 


(PL + Fe Ka) X xabs X xabs 
Xncutral N H (Model E) 


1.90±0.02 


8.37±0.20 


0.61±0.02 
1.27±0.20 


0.36±0.04 
2.33±0.05 


0.05±0.01 


6.41±0.05 




S5 -20 




438/205 


(PL+refl) xxabsxxabs 
Xncutral N H (Model F) 


1.93±0.02 


8.60±0.30 


0.61±0.02 
1.24±0.20 


0.3B±0.04 
2.31±0.06 


0.05±0.01 








0.39±0.10 


443/207 



Table 3. Parameters of fits to the RGS data, compared with EPIC model E. Errors are at 90 % confidence (Ax 2 
interesting parameter. 



2.71) for one 



Instrument 


Warm Nh 


log? 


Neutral N H 


X 2 /dof 




(10 22 cm~ 2 ) 




(10 22 cm~ 2 ) 




EPIC 


0.61±0.02, 1.27±0.20 


0.36±0.04, 2.33±0.05 


0.05±0.01 


438/205 


RGS 


0.20±0.10, 1.31±0.20 


0.50±0.30, 1.62±0.10 


0.12±0.02 


246/231 



that the region that produces the broadest Ha component 
could also be the location of the Fe Ka line emitting gas, if 
this latter emission is due to a single line. 

To start to build up a picture of the AGN as a whole, 
we use the relation in Krolik & Kriss (2001) to estimate the 
approximate distance to the inner edge of the torus due to 
dust sublimation and photoionisation, Rtorus, in cm, from 
the continuum source with ionising luminosity L, in erg s _1 : 



-^torus, cm ^ 3 x 10" 



(1) 



Using the 1-1000 Rydberg luminosity from the EPIC 
best-fit model, i.e. L = 4.3 x 10 44 erg s _1 , we obtain 6.2 x 10 18 
cm (~ 2 pc) for the inner radius of the torus. 

To place a distance on the Fe Ka line formation region 
from the black hole, we use a black hole mass of 3xl0 7 
Mq. This value is intermediate to the Schwarzchild and Kerr 
black hole masses measured for H 0557— 385 by Rokaki & 
Boisson (1999), who assume that H/3 emission lines come 
from the disc. We use the FWHM of 16,000 km s _1 of the 
Fc Ka line as v in the equation: 



GM 
~R~ 



(2) 



to obtain the distance, R, of the line formation region from 
the black hole, assuming the line emitting gas moves in a 
Keplarian orbit: we obtain a distance of 2xl0 15 cm. The 
Schwarzchild radius for the black hole, 2GM/c 2 , will be 
9xl0 12 cm. If we take the inner edge of the accretion disc 
to be 6GM/c 2 , then for H 0557-385 this is 3xl0 13 cm. For 
the Ha line we find distances of 4xl0 17 cm, 5xl0 16 cm, 



and 3xl0 15 cm for the narrowest, intermediate and broad- 
est components respectively. 

4.1 Where did the warm absorbers originate? 

We now try to determine the distances of the warm ab- 
sorbers from the continuum source, and thus find some clue 
as to their origins. 

We assume the warm absorbers to be a continous, con- 
stant velocity wind. We use R; as the launching radius of 
the wind; for such a wind, N H = niRif and £ = L/niR 2 . 
L is the 1-1000 Rydberg luminosity (4.3xl0 44 erg s -1 from 
the EPIC best-fit model) , n; is the density at Ri and / is the 
volume filling factor of the gas. Substituting to eliminate ni , 
we obtain 
r ZNnRt 



f 



(3) 



Taking the upper limit of the volume filling factor, / = 
1, and including the 90 per cent confidence limits on the 
other parameters, we obtain upper limits of Ri < 1 x 10 21 
cm to the higher-^ absorber, and Ri < 3xl0 23 cm to the 
lower-£ absorber. Unfortunately, the locations of the warm 
absorbers cannot be constrained any better from the data 
available. 



4.2 A torus origin for warm absorbers? 

From the relation for the ionisation parameter, £ = L/nr 2 , 
we infer that for constant density and £, r scales with L 1 / 2 . 



XMM-Newton observations of the Seyfert 1 AGN H 0557-385 5 




0.5 1 2 5 10 0.5 1 2 5 10 

Energy (keV) Energy (keV) 



1 _ 1 1 1 






c 


: % = 637/209 




: - ^ ... 'J 


i 





Energy (keV) 



D 



X = 487/208 



"*V , i '^/* ' i | ' | I ''>''' i /' | V* 



2 

Energy (keV) 





0.5 1 2 5 10 0.5 1 2 

Energy (keV) Energy (keV) 



Figure 1. Spectral fits of EPIC data, plotted in the observed frame on logarithmic scales. The plots below each fit show the data/model 
ratio. A: Power law, fitted over 2.3-10 keV, plotted back to 0.3 keV. B: Power law and one warm absorber. C: Power law and two warm 
absorbers. D: Power law, two warm absorbers and neutral absorption. E: as D but including an iron Ka line. F: as D but with a reflection 
component. 



Equation expresses the same scaling relation, for the in- 
ner edge of the torus. Here, r is the distance of the ionising 
source from the absorbing gas. From this we can infer that 
if warm absorbers originate at the inner edge of the torus 
then they should possess similar values of the product of ion- 
isation parameter and density, independent of their ionising 
luminosities. 

We re-write Equation Q as fltorus ~ KyL with K 
= 3xl(T 4 . Then re-arranging Equation [3] equating Ri to 
fltorus, and substituting for i?torus, we obtain 



L^i^ (4) 

Taking logs of both sides then working through the al- 
gebra, we obtain 

Jo 5 L~2 1og£iVH-2 1og-^ (5) 

This implies a linear relationship between log L and 
log I; Nh with a gradient of 2 and an intercept of —2 log -L. 
A list of logJVjj, log£ and logL values for warm absorbers 
observed in 14 AGN is given in Tables 2 and 4 of B05. These 
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Figure 3. The RGS spectrum, plotted in the observed frame with the best-fit RGS model. Some of the main absorbing ions that belong 
to each phase are labelled. The data/model ratio is also shown. 



authors assume a constant velocity wind in considering the 
energetics of warm absorbers. In Fig. we plot log L versus 
log £Nn for the AGN in B05 along with the parameters we 
derive for H 0557-385 from the RGS data. 

We can use Equation to calculate the filling factor 
/ when all the other terms are known. The values of L, 
A^h and £ have been determined from our observations of 
H 0557-385 and are listed in B05 for 14 other AGN. The 
filling factors for each absorber phase are listed in Table [3] 



The absorber of PG 0844+349 and one of the absorbers 
in PG 1211+143 would have filling factors larger than 1 if 
they originate at the inner edge of the torus, so they cannot 
originate as torus winds. This would also hold if we consider 
accelerating rather than constant velocity winds, as in this 
case Njj < niRif and the apparent value of / would be 
larger than that indicated in Table [I] Pounds et al (2003a, 
2003b) and B05 argue that the absorbers in these objects 
originate as accretion disk winds; from Fig. |S] these must 
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20 

Wavelength A 

Figure 4. The RGS best-fit model, plotted in the observed frame, for the two absorber phases. Black: high-£ phase. Grey: low-£ phase. 
The top dotted curve is the unabsorbed continuum; the lower curve, shown in bold black, includes the neutral gas. The Galactic column 
is not included in any of the models in this plot. 



20 




28 



Figure 5. Log values of ionising luminosity, versus the log of the product of ionisation parameter and column density, for 14 AGN. 
H 0557—385 is shown as starred points, along with the errors. Data for all other AGN are from B05. The warm absorbers in PG 0844+349 
and PG 1211+143, believed to be accretion disk winds, are shown as triangles. Some different filling factors are labelled. 



originate closer than the torus, so the results are consistent 
with an accretion disk wind. 

From Fig. E|and Table[I]there are 6 absorbers for which 
we compute very low filling factors (/ < 1CP 4 ). These ab- 
sorbers have the smallest values of log £ in the sample, rang- 
ing from to 0.68. The AGN in question are H 0557-385, 
MR 2251-178, IRAS 13349+2438, NGC 4593, NGC 5548 
and Ark 564. 

In Fig.EJwe plot the filling factors for the AGN as listed 
in Table |1J versus those found by B05. The advantage of 
our method for calculating filling factors compared to that 



of B05 is that our method only requires the logA^jj, log£ 
and log L values for warm absorbers, and no dynamical in- 
formation is required. In contrast, the method of B05 used 
the outflow velocities of the warm absorbers, and involves 
calculating the momentum absorbed and scattered by the 
warm absorber. So two different methods exist for finding 
the filling factor of the warm absorber. 

Unsurprisingly, Equation |H| assigns much larger filling 
factors to the accretion disk wind absorbers than B05. For 
the low-£ absorbers, the two methods produce inconsistent 
results. For the remaining objects, the absorbers with log 
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£ > 0.7 are located in a different part of Fig. [(J to the 5 
lower-£ absorbers, which are found on the left-hand side of 
the plot. The absorbers with log £ > 0.7 have -4 < log / < 
in both methods. Their mean log / from B05 is -1.8±0.2 
with a standard deviation of 0.7, and their mean log / from 
this work is -2.1±0.2 with a standard deviation of 0.8 where 
the errors are 1 a. Therefore the two methods yield consis- 
tent values for (log /} . It is worth comparing the standard 
devations of these results with the uncertainties inherent in 
the two methods. For our method we estimate the overall 
uncertainty on / to be approximately 0.5 dex (a sum in 
quadrature of 0.15 dex from L ion , 0.2 dex from f, 0.3 dex 
from Nh and 0.3 dex from K, the latter from the uncer- 
tainty in the distance of the inner edge of the torus). For 
the method of B05, we estimate an uncertainty on / of ap- 
proximately 0.6 dex (a sum in quadrature of 0.3 dex from 
Li 0n , 0.2 dex from £ and 0.5 dex from the absorbed and 
scattered momentum). Hence, the standard deviations on / 
derived from the two methods are comparable with the er- 
rors we estimate for the two methods. For the absorbers with 
log £ < 0.7, the mean log / from our method is -5.1±0.1 and 
that of B05 is -2.6±0.3. 

If we relax our assumption that the absorbers originate 
at the inner edge of the torus, we can determine the values 
of Ri that would bring the two methods into agreement, by 
substituting the filling factors of B05 into Equation[3] These 
values of Ri are given in column 8 of Table 01 Taking the 
standard deviations on / found from the high-£ group of 
absorbers as the uncertainties inherent in the two methods, 
and summing them in quadrature, we obtain an uncertainty 
of ~ 1 dex. This explains why some of the values of Ri deter- 
mined are very large (the low-£ absorbers in MR 2251-178 
and Ark 564 are estimated to be greater than 3 Kpc distant 
from the nucleus). The mean of log Ri for the low-£ ab- 
sorbers is 2.63, so a typical distance for these absorbers is ~ 
400 pc, which is similar to the scale of the narrow-line region 
(Osterbrock 1993, Capetti et al 1999). At a similar distance 
from the continuum source, a photoionised cloud is seen in 
emission in X-ray observations of the prototypical Seyfert 2 
galaxy NGC 1068 (Brinkman et al 2002, Ogle et al 2003). 
The low-£ absorbers could be equivalent to this photoionised 
cloud in NGC 1068, but seen in absorption rather than in 
emission. It is unclear whether such absorbers originate from 
the torus or from other, more distant, material. 

4.3 Where is the neutral gas component? 

We find that the RGS data of H 0557— 385 require absorp- 
tion from neutral gas local to the galaxy, of column density 
1.2x 10 21 cm -2 . But where is this gas in the AGN structure? 
Fairall, McHardy & Pye (1982) and Rafanelli (1985) found 
the optical emission lines of H 0557— 385 heavily reddened. 
Therefore if we assume that dust is mixed with the neutral 
absorber, then both must lie outside the broad line region. 
There are therefore three main possibilities for the location 
of the neutral gas: the AGN itself, the disk of the host galaxy, 
or somewhere in the host galaxy but outside the disk. 

Scenario 1: The first possibility is gas within the AGN 
itself. If the warm absorber is a wind formed by evaporation 
from the dusty torus, then the neutral gas could be from 
the torus. Ogle (2003) suggests that the torus has a role 
in defining the shape of the ionisation cone. As we must be 



viewing the continuum source through the ionisation cone to 
see a warm absorber, we cannot be looking directly through 
the side of the torus. However, it could be that we are seeing 
cold absorption through a small cross-section of the torus if 
our line of sight grazes the edge. 

Scenario 2: The second possibility is the disk of the 
host galaxy. H 0557—385 is classified as SOa in the Nasa 
Extragalactic Database, between a spiral and a lenticular 
galaxy. It is observed at an inclination angle i = 75° from 
the normal (Crenshaw & Kraemer 2001). If the host galaxy 
of H 0557—385 has similar column density values to other Sa 
galaxies, then the neutral absorption we observe could orig- 
inate in the host galaxy of the AGN. The average neutral 
hydrogen column density for an Sa galaxy viewed face-on, 
N H± , is 1.4xl0 21 cm~ 2 (Broeils & Woerden 1994). Using 
the relation iVn = iVux/cos i, we obtain an expected value 
of 5.3xl0 21 cm -2 for Nr. This is larger than the neutral 
column we observe in H 0557— 385, so it is plausible the neu- 
tral gas in H 0557—385 could originate in the disk of the host 
galaxy. 

Scenario 3: The third possibility for the origin of this 
neutral gas is outside the AGN and above the disk of the 
galaxy, but within the structure of the galaxy. It could be in 
the form of a dust lane. From HST images, Malkan, Gorjian 
& Tarn (1998) find that there are dust lanes at distances 
of hundreds of parsecs in Seyfert 1 and 2 AGN. These dust 
lanes have little or no connection with the central engine and 
could have column densities up to ~ 10 23 cm -2 (Matt et al 
2000). Such a dust lane could therefore plausibly supply a 
column of ~ 10 21 cm -2 . 

Given these three scenarios, we consider that all possi- 
bilities are plausible. 

Since some of the values of R; (the launching radius of 
the wind) for different AGN in Table H are on similar scales 
to those of Scenario 3, it may be that the cold gas can exist 
on similar scales to the warm absorbers, this is plausible 
given the small filling factors of the warm absorbers. 

4.4 Conclusions 

From our XMM-Newton RGS observations of the Seyfert 1 
AGN H 0557— 385 we find that warm absorption is present 
in two phases, a low-^ phase with log £ of 0.50 and a col- 
umn of ~ 10 21 cm -2 and a high-£ phase with log ionisation 
parameter £ of 1.62 and a column of ~ 10 22 cm -2 . We de- 
tect an unresolved transition array of (M-shell) iron in both 
phases of the absorber. 

We infer a picture of the different scales of these com- 
ponents. The locations of the warm absorbers are not well 
constrained ($: 10 21 - 3x 10 23 cm from the central continuum 
source), compared to the molecular torus distance of ~ 10 18 
cm from the source of the ionising luminosity. The Fe Ka 
line and a broad Ha line are both likely to originate at ~ 
10 cm from the black hole. 

The neutral gas observed could originate from the AGN 
torus, the host galaxy disk, or a lane or cloud of gas above 
the disk. 

We establish a new method of determining the volume 
filling factors in warm absorbers, assuming that the ionised 
absorbers originate as an outflow from the torus. We derive 
the volume filling factors for the Blustin et al (2005) sam- 
ple of AGN and compare them with those published. We 
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Figure 6. The filling factors from this paper plotted versus the filling factors derived from B05, which are translated from percentages 
into a fractional scale. The warm absorbers in PG 0844+349 and PG 1211+143, believed to be accretion disk winds, are shown as 
triangles. A line indicating equivalence between our results and those of B05 is shown. 



find that they have small volume filling factors 1%). We 
find reasonable agreement with the method of Blustin et al 
(2005). However, a group of five absorbers have filling fac- 
tors derived in this paper that are much smaller than those 
found by B05. These absorbers all have values of log £ that 
are < 0.7, and appear to lie out in the narrow-line region, 
typically hundreds of parsecs. 
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Table 4. Comparison of volume filling factors (f) between this work, and B05. The B05 filling factors are translated from percentages 
into a fractional scale, 'w.a.' stands for warm absorber. Rtorus is the distance of the inner torus edge, and i?; is the launch radius 
calculated for the warm absorbers using the filling factors from B05. The uncertainty on R; is approximately ±1 dex, see Section. 14.21 
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